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AN  EQUATION  OF  STATE  FOR  SODIUM^ 


P.(  Browning 


A  revised  equation  of  state  for  sodium  is  presented.  Changes  in  reference  data 
and  improvements  in  computational  technique  to  those  previously  used  are  described. 
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\  Introduction 

m- 


( EOS ) ^ 


The  equation  of  state  (EOS)  for  sodium  which  has  been  enployed  in  assessments  of 
hypothetical  accidents  in  liquid  metal  cooled  fast  breeder  nuclear  reactors  has  been  in 
use  for  some  years  in  the  British  programme.  During  this  time  some  important  experi¬ 
mental  reference  data,  upon  which  the  EOS  is  based,  have  been  revised.  The  purpose  of 
this  report  is  primarily  to  update  the  sodium  EOS  by  incorporating  these  revised  data. 

In  addition,  a  number  of  improvements  have  been  made  in  the  calculational  technique  used 
in  deriving  properties  in  the  single  phase.  These  refinements ,  which  have  indicated 
numerical  errors  in  the  earlier  EOS  output  ,  have  improved  the  precision  of  the  reported 


he  Equation  of  State 


A  brief  resume 


As  in  Alderson's  original  report,  the  P-V  surface  (Fig.  1)  is  divided  into  four 
regions: 

(i)  The  compressed  or  subcooled  liquid 

(ii)  The  expanded  or  superheated  gas 

(iii)  The  two  phase  coexistence  region 

(iv)  The  supercritical  region. 

Evaluation  of  thermodynamic  properties  within  each  of  these  regions  is  based  upon 
a  quantification  of  the  envelope  denoting  the  boundary  between  the  two-phase  coexistence 
region  and  the  single  phase  regions. 

The  equilibrium  envelope  is  effectively  characterised  once  all  the  equilibrium 
properties,  shown  in  Table  I,  are  defined  as  a  function  of  temperature.  With  respect  to 
the  determination  of  the  position  of  the  equilibrium  envelope,  the  thermodynamic  proper¬ 
ties  within  the  coexistence  region  are  evaluated  by  application  of  the  lever  rule. 

The  properties  in  the  compressed  liquid  and  the  expanded  gas  regions  are  evaluated 
on  the  basis  that  the  thermal  pressure  coefficient,  Tv  -(§?)„«'"  Appendix  1 ) ,  is  known 

throughout  these  regions.  As  no  data  are  available  on  the  behaviour  of  xv  in  the  single 
phase,  the  assumption  employed  in  both  this  and  earlier  work  is  that  Yy  is  only  a  function 
of  volume.  Consequently,  provided  it  can  be  evaluated  at  one  volume,  that  value  can  then 
be  applied  at  the  same  volume  for  any  other  temperature. 

In  general,  the  evaluation  of  data  for  the  subcritical  single  phase  region  involves 
integration  of  the  relevant  quantity  between  the  point  of  interest  and  the  appropriate 
point  on  the  saturation  envelope.  Since,  above  the  critical  temperature,  the  saturation 
envelope  is  undefined,  a  further  assumption  needs  to  be  employed  in  order  to  characterise 
the  supercritical  region.  The  assumption  is  that,  at  supercritical  temperatures  (TSUp)> 
the  internal  energy  on  the  critical  isochore  V,,  is  directly  proportional  to  the  tempera¬ 
ture  difference  (Tg  -  Tc).  The  constant  of  proportionality  is  based  on  an  averaging 
of  the  internal  energies  on  the  critical  isochore  for  some  small  temperature  range 


directly  below  T  and  therefore  generating  an  effective  specific  heat  C  (V  ).* 

C  VC. 

In  the  absence  of  any  supercritical  data  for  sodium  this  treatment  does  at  least 
allow  the  supercritical  region  to  be  quantified,  but  the  method  cannot  be  rigorously 
justified.  This  second  assumption  effectively  establishes  the  critical  isochore  as  a 
reference  curve  in  the  same  sense  as  the  coexistence  envelope.  A  future  development 
which  would  obviate  the  use  of  a  critical  isochore  reference  line  would  be  to  establish 
the  position  of  a  shock  Hugoniot  within  the  compressed  liquid.  This  would  provide  an 
experimental  reference  line  from  which  extrapolations  could  be  made  with  greater 
confidence. 

The  above  description  summarises  the  approach  which  has  been  used  to  obtain  an  EOS 
code  used  for  sodium.  The  techniques  used  are  exactly  those  employed  by  Alderson^. 

As  the  new  EOS  is  generated  in  exactly  the  same  wav  there  is  no  detailed  description  of 
the  determination  of  the  reference  curves  or  of  the  evaluation  of  the  thermodynamic 
properties,  but  simply  the  equations  used  are  presented.  It  has  been  left  to  the  reader 
to  consult  the  original  work  of  Alderson  for  a  full  description.  The  sections  which 
follow  describe  the  revisions  which  have  been  made  to  that  original  treatment. 


3.  The  Revised  Critical  Data 

Accurate  knowledge  of  the  critical  temperature,  volume  and  pressure  is  crucial  to  a 

(2) 

confident  construction  of  the  equilibrium  envelope.  The  measurements  of  Bhisa  have, 
for  some  time,  been  accepted  as  the  most  reliable^  values  for  the  critical  constants 
of  sodium.  These  are  shown  in  Table  II,  along  with  those  previously  used  in  Alderson's 
equation  of  state. 


The  adoption  of  a  new  set  of  critical  constants  alters  many  of  the  equations 
describing  the  equilibrium  properties  on  the  coexistence  envelope.  In  this  treatment 
of  the  sodium  EOS,  Alderson's  original  equations  are  used  up  to  1640K.**  Above  this 
terrqierature  equations  of  exactly  the  same  form  are  used  except  that,  in  order  to  accommo¬ 
date  the  new  critical  constants,  the  coefficients  of  the  equations  necessarily  change. 


In  general,  the  coefficients  of  an  equation  representing  a  given  property  above 
1640K  can  be  determined  with  a  knowledge  of  the  critical  point  values  themselves,  and  by 
the  continuity  constraints  at  1640K,  e.g.  the  three  coefficients  of  the  expression  for 
the  vapour  pressure,  P  : 

Sal 


P  „  (T  >  1640K) 
sat 


(1) 


are  fixed  by  a  knowledge  of  the  critical  pressure,  the  pressure  at  1040K  and  the  slope  of 
the  vapour  pressure  with  respect  to  temperature  at  1640K,  these  last  two  values  being 
obtained  from  the  vapour  pressure  equations  at  lower  temperatures.  Occasionally  a 
continuity  constraint  has  no  thermodynamic  significance,  but  is  purely  a  mathematical 


*  The  variation  of  internal  energy  with  t  em[>erat,uro  along  Vc  and  below  Tc  is  determined 
by  application  of  the  lever  rule  to  the  internal  energies_of  liquid  and  vapour. 

From  the  resulting  internal  energies,  an  'average'  slope  CylV,,)  can  be  established, 
which  is  then  used  at  temperatures  above  Tc. 

**  As  distinct  frcm  I644.28K  which  has  no  significance  apart  from  being  equivalent  to 
2 500°F . 
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convenience.  In  the  case  of  the  vapour  pressure,  however,  the  slope  dP/dT  is  related 
to  the  heat  of  vaporisation.  This  must  vary  smoothly  across  the  l<>40K  boundary  as  there 
is  no  phase  change  for  condensed  sodium  at  this  temperature. 

The  revised  equations,  which  replace  those  of  Alderson  in  the  new  code,  are  given 
in  Appendix  2. 

4.  Improvements  in  the  Precision  of  Calculations 

An  improvement  in  the  accuracy  of  characterising  the  equilibrium  envelope  must 
result  from  an  improvement  in  an  estimate  of  the  critical  constants.  The  precision  with 
which  the  thermodynamic  properties  are  determined  away  from  this  reference  envelope  is, 
however,  also  dependent,  to  various  extents,  upon  the  mathematical  techniques  involved 
in  their  evaluation.  In  applying  the  equations  of  Appendix  2  to  the  determination  of 
single  phase  properties  it  becomes  apparent  that  both  nunerical  integration  and  numerical 
differentiation  are  unavoidable. 

The  eventual  formulation  of  the  equations  of  Appendix  2  into  an  EDS  for  sodiun  was 
achieved  independently  from  the  code  used  by  Alderson.  In  this  way  the  new  EOS  provided 
a  useful  test  of  the  precision  of  the  Alderson  code. 

Comparison  of  the  results  from  the  two  codes  revealed  differences  in  the  prediction 
of  certain  parameters  at  low  temperatures  (T  <  700K)  which,  by  virtue  of  being  described 
by  the  same  equations,  should  not  have  existed.  Investigation  of  both  codes  revealed 
that  the  numerical  differentiation  technique  used  by  Alderson  was  inadequate  at  low 
temperatures.  This  is  discussed  fully  in  Appendix  3;  at  this  point  it  is  sufficient  to 
say  that  the  values  of  and  used  by  Alderson  for  the  equilibrium  vapour  below  700K 
are  approximately  twice  those  given  in  this  EOS.  This  has  important  consequences,  as 
will  be  seen  in  Section  5. 

A  further  comparison  of  the  two  codes  showed  a  disagreement  in  the  high  specific 
volume  values  of  internal  energy  (U)  and  enthalpy  (H)  along  isotherms  near  and  above  Tc< 
This  is  again  discussed  in  Appendix  3  but  essentially  the  rapid  change  of  U  and  H  at  high 
volume  is  a  consequence  of  an  inefficient  numerical  integration  technique  used  by 
Alderson.  The  internal  energies  and  enthalpies  ought,  in  fact,  to  tend  to  their 
equivalent  perfect  gas  limits  at  high  volumes  for  a  given  temperature,  not  diverge  as 
shown  in  Aiders  on's  tables. 


5.  Results  and  Discussion 

The  equations  employed  in  this  EOS  are  given  in  Ap|>endix  2.  The  EOS  itself  is  given 
in  microfiche  form  contained  within  this  report.  The  structure  of  the  EOS  is  in  the  same 
format  as  the  original**^.  Appendix  4  lists  the  data  contained  within  the  microfiche. 

Some  comment  is  required  on  whether  or  not  this  reassessed  EOS  for  sodium  represents 
an  improvement  on  Alderson' s  EOS  both  from  a  purely  thermodynamic  consideration  and  for 
the  analysis  of  hypothetical  accidents.  A  more  accurate  set  of  critical  constants  must 
improve  the  confidence  placed  upon  the  description  of  sodium  liquid/vapour  equi libriun. 

The  low  temperature  specific  heats  of  the  vapour  are  greatly  improved  by  the  implementation 
of  a  more  efficient  calculat  ional  teclinique.  In  this  sense,  when  equilibrium  data  are 
required  in  hypothetical  accident  analysis  the  results  of  such  an  analysis  must  be  improved. 


,  *#**■.  .  *  .x* in 


When  single  phase  data  are  required  for  analysis  of  hypothetical  accidents  the  concern 
over  the  precision  of  the  extrapoiative  technique  is  the  same  as  for  the  former  EOS. 

Some  attempt  to  quantify  the  precision  of  such  predictions  can,  however,  be  made, 
depending  upon  the  extent  into  the  single  phase  regions  to  which  the  extrapolations  are 
applied.  Calculations  in  the  low  temperature  expanded  gas  region  are  obviously  quite 
adequate.  Below  1400K  the  vapour  is  virtually  ideal.*  For  an  ideal  gas  the  Yy 
extrapolation  is  no  longer  an  approximation,  it  is  exact.  Consequently,  below  1400K 
the  expanded  gas  extrapolations  are  as  accurate  as  the  equilibrium  vapour  reference  line 
upon  which  they  are  based. 

Above  1400K,  up  to  ca.  I800K,  the  "ideality  factor"  is  still  greater  than  ca. 

0.7.  As  is  later  demonstrated  in  Appendix  1,  Yv  for  the  vapour  is  evaluated  from  a  van 
der  Waals  type  treatment.  At  ca.  1800K  the  specific  volume  of  the  vapour  is  still 
sufficiently  large  to  allow  a  van  der  Waals  type  description  of  the  onset  of  non-ideality. 
Consequently,  the  extrapolations  into  the  expanded  gas  between  14O0K  and  1800K  are  a  fair 
approximation  to  reality. 

The  use  of  the  Yv  extrapolation  in  the  compressed  liquid  region  does  not  benefit 

from  a  limiting  comparison,  as  does  the  expanded  gas  •-  )ierfect  gas.  The  temperature 

(4) 

insensitivity  of  y  (constant  along  an  isochore)  is  an  observed  feature  for  some 
materials.  A  strong  indication  that  the  -ame  constraint  is  valid  for  sodium  is  demon¬ 
strated  by  comparing  the  Yv  extrapolation  with  one  based  on  the  so-called  Gruneisen 

v  /3p\ 

extrapolation.  This  asserts  that  the  Gruneisen  parameter  =  (  l/V)  V  is  a^so 

constant  along  an  isochore.  If  both  extrapolations  give  rise  to  the  same,  or  very 
similar,  extrapolated  results,  our  confidence  will  increase.  Furthermore,  if  the 
equilibrium  Cy  (for  both  phases)  can  be  shown  to  he  reasonably  insensitive  to  temperature 
(see  Appendix  5),  then  the  extrapolation  must  be  considered  reasonable.  Alderson 
demonstrated  this  to  be  the  case  for  sodium  except  for  the  low  temperature  vapour.  As 
was  mentioned  in  Section  4,  Alderson's  values  for  Cv  for  his  low  temperature  vapour  were 
over-estimated.  Fig.  3  shows  that  the  Cy  agreement  is  improved  over  that  of  the  previous 
work . 

The  values  of  Cy(liq)  and  Cy(vap),  whilst  being  ca.  \5%  different,  do  exhibit  the 
same  trends  over  the  range  600  <  T  <  2300K.  Fig.  5  gives  an  indication  of  the  agreement 
between  the  Yv  and  Ty  extrapolations.  From  this  we  can  infer  that  up  to  2300K  the 
techniques  predict  virtually  the  same  results. 

For  2300  '  TOO  <  Tc,  the  variation  of  Yy  with  temperature  no  longer  benefits  from 
comparison  with  either  a  perfect  gas  limit,  or  a  constancy  in  C  .  Justification  for 
continuing  with  a  constant  Yy  at  these  temperatures  rests  purely  with  observation  of  the 
invariance  of  Yy  with  temperature  for  other  materials.  Above  Tc  the  extrapolations,  as 
mentioned  in  Section  2,  are  based  upon  the  assumption  of  a  constant  Cy**  on  the  critical 


*  PV/RT  ?  0.85  for  T  «  HOOK. 


**  Cy  is  not  a  constant  on  an  isochore  for  subrritical  temperatures. 
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isochore  above  T  (not  exactlv  at  T  as  C  =  no) .  This  constraint  is  perhaps  the  least 
easily  defended  in  the  quantification  of  this  type  of  EOS.  The  constancy  of  Cy  along  the 
critical  isochore  is  consistent  with  a  Yv  description.  However,  Tc  is  a  point  of 
singularity  with  respect  of  Cv.  The  adoption  of  a  Cy(T  >  T  )  having  a  value  equal  to 
the  average  value  obtained  by  applying  the  lever  rule  on  the  critical  isochore  below  Tc 
is  an  intuitive  assignment  in  the  absence  of  any  experimental  guidance,  and  therefore 
represents  a  best  approximation  at  this  time. 

In  conclusion,  the  reported  EOS  is  considered  an  improvement  over  Alderson's 
original  tables  and  its  adoption  will  improve  the  thermodynamic  description  of  sodium. 
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TABLE  I 


Equilibrium  Properties  Needed  to  Characterise  Fully 
the  Saturation  Envelope 


The  following  equilibrium  quantities,  as  functions  of  temperature  up  to  T 
are  sufficient  to  characterise  fully  the  EOS  of  sodiun.  ( 

Vapour  pressure 

Liquid  and  vapour  densities 

Liquid  and  vapour  enthalpies 

Liquid  and  vapour  adiabatic  compressibilities 

Vapour  thermal  pressure  coefficient 

Those  quantities  not  mentioned  in  this  list  are  directly  calculable  from 
the  above  quantities. 


TABLE  II 


The  Old  and  Revised  Critical  Constants  for  Sodium 


APPENDIX  1 


The  Evaluation  of  Single  Phase  Properties 
The  thermal  pressure  coefficient,  Yv ,  is  defined  as 


From  the  considerations  of  Section  2  the  pressure  in  the  single  phase  (sp)  region  is 
directly  obtained  from 


sp 


=  P  „  +  Y  (T  -  T  . ) 
sat  v  sp  sat 


(1.1) 


where  Yv  is  evaluated  at  T^at  on  the  two  phase  envelope.  As  is  discussed  later  in  this 
appendix,  the  behaviour  of  the  thermal  pressure  coefficient  determines  all  the  single  phase 
thermodynamic  properties.  The  calculation  of  Yv  is  therefore  crucial. 

Evaluation  of  Yv  along  the  liquid  saturation  line  is  described  in  Reference  I.  The 
approximation  that  Yv  is  only  a  function  of  volume,  thereby  allowing  extrapolative  calcula¬ 
tions  to  be  performed,  is  based  on  observation.  No  further  comment  on  the  extrapolation 
of  Yv  (sat.liq.)  will  be  made  here.  The  extrapolation  of  Yy( sat.vap.)  does  benefit  from 
comparison  of  Yy  for  both  a  perfect  gas  and  a  slightly  "imperfect"  gas.  The  expression 
used  to  evaluate  Yv  (sat.vap.)  in  this  work  is  of  the  form: 

Yy(  sat.  vap. )  =  - - - - — g  (1.2) 

sat.vap. 


where  A  and  B  are  constant. 

If  we  consider  the  perfect  gas  equation 

PV  =  RT 


then: 


(1.3) 

(1.4) 


Comparing  (1.2)  and  (1.4)  we  see  that  for  high  volumes  the  general  form  of  (1.2)  tends 
towards  (1.4)  (i.e,  B  can  eventually  be  neglected).  Titus  for  (1.2)  to  be  a  reasonable 
equation  to  use  the  coefficient  A  must  be  close  to  the  gas  constant,  R.  The  value  is 
shown  in  Appendix  2  and  is  indeed  seen  to  be  close  to  R.  Thus  the  limiting  form  of  (1.2) 
is  correct.  As  V  decreases  then  B  can  no  longer  be  ignored.  The  origin  of  B  can  be 
demonstrated  using  van  der  Waals'  equation: 

("•■?)•  (v " 3)  ‘  RT  "•5) 
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for  small  departures  from  ideality  (1.5)  reduces  to 


P(V  -  P)  =  RT 


/ap\  _  r 

Yv  ~  l  0T  )  =  V  -  P 

'  ' v 


Comparing  (1.6)  with  (1.2)  we  again  see  that  A  can  be  associated  with  R,  and  B  with  3,  in 
the  van  der  Waals  equation. 

The  van  der  Waals  equation  is  most  useful  as  a  qualitative  description  of  the  onset 
of  non-ideality  in  a  gas.  In  practice  it  has  little  quantitative  application  except  at 
very  high  volumes.  The  choice  of  this  type  of  expression  to  describe  Yv  (sat.vap.), 
right  up  to  Tc,  lies  with  Alderson^  .  The  author  was  aware  of  the  inherent  inadequacies 
and  chose  to  allow  both  A  and  B  (eqn.  (1.2))  to  vary  from  their  high  volume  values,  to 
afford  a  better  description  near  Tc,  rather  than  use  a  more  appropriate  expression. 

For  the  purpose  of  this  EOS,  evaluation  of  single  phase  properties  reduces  to  that  of 
obtaining  values  for  the  internal  energy  and  entropy.  As  mentioned  in  Section  2,  this 
involves  integration  of  the  appropriate  function  with  respect  to  a  point  on  the  reference 
curve  (Tq,Vo).  The  energy  in  the  condensed  liquid  region  is  evaluated  using 


U(To,v,)  =  U<To,V( 


V 


(Ty  -  P)  dV 
'v  'sat 


Equation  (1.7)  is  best  illustrated  by  reference  to  Fig.  3.  The  integration  in  (1.7)  is 
evaluated  numerically  by  evaluating  (TY  -  P)  o4.  at  an  adequate  number  of  points  along  the 

V  Sat- 

two  phase  envelope  between,  and  including,  VQ  and  Vj.  The  entropy  at  the  same  point  is 
evaluated  using 


S(T0’Vi ) 


S(To’Vo> 


(Y  )  .  dV 

'v  sat 


where  the  quantities  in  (1.8)  are  entropy  equivalents  of  the  energies  in  (1.7). 

The  evaluation  of  these  quantities  in  the  expanded  gas  is  achieved  in  exactly  the  same 
fashion.  Up  to  the  critical  isotherm,  the  energy  and  entropy  are  evaluated  in  the  above 
manner.  Once  into  the  supercritical  region  they  are  extended  using  Cv(Vc),  introduced  in 
section  2.  Fig.  4  illustrates  the  technique  for  evaluating  supercritical  energies; 
entropies  are  evaluated  in  an  analogous  manner,  the  entropy  along  the  critical  isochore 
being  evaluated  bv 

T 

S(To’Vc)  =  S(Tc’Vc)  +  ev(Vc>  *"  (f)  (,-9) 


and  the  energy  bv 


U(To’V  =  U(TC’V  +  5v(VcMTo  ”  V 


(1.10) 


-10- 


I 


\ITL\UIX  2 


Cpdal  < ■< I  Ennui  inns  for  Na  for  l_  I  h)40K 


Critical  Constants: 


Saturated  Vapour  Pressure: 


25.98  MPa 


V  4.885  /  10 

c 


P  (T)  =  Ae  '  T  (MPa) 


A  =  48(1 .714 


I  1419.(18 


-  -0.2085587 


Saturated  Liquid  Volume: 


V  (T)  = 


Ini’  \n) 


1  *  A(  1  -  HT) 


5.995205755  '  10 


0.489129 


Saturated  Vapour  Volume: 


VSV(T)  -  Vsl(T)  +  (P,.  (0.2085587  .IHH.wfl  ^ 


AH  is  the  onthnlpv  of  vaporisation  anil  i-  xiven  lielow. 
vap 


Saturated  Linuid  Entropy: 


S  .jIT)  =  S^(  K140K) 


7  C  , ( T)  dT 


(k.l  kK.K) 


(T»  ,„o 

C.  (T)  - ^ - V  .  -p-  -  Id  (k.l  kg.K) 

si  \  <n  s i  <n  ; 

v  4' 


dll  .IT)  A.D.CI  I  -  T  I  I 
-  I  _ 

dT  2. ( T  -  T) 


(Ml  k«.K) 


8.51525  '  It) 


1.9  152  '  IO 


-I  I- 


and 


nil  ere 


dP 

Vd? 


si 


Ae  'c.T*-0-1^  +  T 

* 


A  =  480.714 

B  =  11419.68 

C  =  -0.2085587 


Saturated  Liquid  Enthalpy: 


where 


Hsl(T) 


A 

B 


A  +  BT  -  4  AH 

2  vap 


2.4423 

8.5125  x  10-4 


Enthalpy  of  Vaporisation: 


where 


“Vap  =  ATC"'VB 

A  =  1.9152  X  10~3 

B  =  0.33699 


Saturated  Vapour  Thermal  Pressure  Coefficient: 


where 


YV(T) 


A _ 

V  x  i03  -  B 
sv 


A  =  .382268 

B  =  0.02959 


(.Note:  the  gas  content  R  =  0.36  16  37  MPa.  m3/kg.K  ') 
Saturated  Liquid  Adiabatic  Compressibility: 


a  _  .  /  T  -  B  \ 

s  ^C<Tc  -  T)J 

where 

A  =  1.985  x  icf4 

B  =  5252.0 

C  =  3.0751 


Saturated  Liquid  Thermal  Expansion  Coefficient: 


a 


si 


where  V  . 

s  l 


B  \  (MPa.K- 1 ) 


(MJ/kg) 


(MJ/kg) 


(MPa/K) 


(MPa'1) 


is  ,  iven  above,  and 
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where 


dVsl  _  ABC( 1  -  BT)C~ 1 
dT  1  +  A( 1  -  BT)C 


(m  /kg.K) 


A 

B 

C 


3.23112 
3.995205  x  icf 
0.489129 


All  other  thermodynamic  quantities  evaluated  in  this  EOS  are  calculable  given  the 
above  equations.  The  reader  is  referred  to  Reference  1  for  these  additional  equations. 


AITIVUX  3 


Improvement  s  in  Ca  leu  tat  i  onal  Technique 


Tun  part-  of  the  data  given  in  Reference  I  are  in  error.  These  are  the  low  t  empr  ra- 
t  lire  -peri  fir  heat,  C  ,  of  the  saturated  vapour,  anti  the  high  volume  internal  energies  in 
the  -ingle  phase  expanded  gas. 


It  is  necessary,  in  calculating  C  for  both  the  liquid  and  vapour,  to  evaluate  the 
-aturated  expansion  coefficients,  :  . ,  for  hoth  phases 


(3.1) 


where  j  represents  either  liquid,  1,  or  vapour,  v.  For  the  liquid  this  may  be  done 
analytically  using  the  equation  given  in  Appendix  2.  For  the  vapour,  owing  to  the  proce¬ 
dure  involved  in  evaluating  Vgv  it  is  much  easier  to  calculate  a  numerically.  In  fact, 
the  evaluation  of  a  ^  is  also  done  numerically  for  convenience,  no  significant  loss  of 
accuracy  being  incurred.  In  the  evaluation  of  this  numerical  differential  it  is  necessary 
to  choose  a  suitable  step  length  in  T.  Alderson  chose  a  5K  step  length.  This  is  ade¬ 
quate  for  most  calculations,  however,  at  temperatures  helo-  700°C  the  specific  volume  of 
the  vapour  increases  rapidly  with  decreasing  temperature.  A  step  length  of  5K  is  too 
large  for  a  reliable  estimate  of  agv  in  this  region.  The  consequence  of  this  was  to 
over  estimate  by  ca.  1001?  at  f>00K.  In  Section  4  it  was  stated  that  this  had  important 
consequences  regarding  the  general  applicability  of  the  EOS.  It  will  he  seen  in 
Appendix  3  that  one  of  the  arguments  supportin’-  the  adoption  of  a  Yy  which  is  only  a 
function  of  volume  is  the  'invariance'  of  Cv  with  temperature,  for  both  liquid  and  vapour. 
In  his  original  work  Alderson  concluded  that  "C  is  approximately  constant  up  to  about 
2500K  apart  from  the  low  temperature  vapour."  Fig.  2  shows  Cv  for  hoth  phases,  evaluated 
from  this  work  and  the  agreement  at  low  temperatures  is  improved,  compared  with  the 
original  work.  in  consequence,  confidence  in  the  Yy  extrapolation  technique  is  increased. 

The  second  area  of  Alderson 's  tabulation  in  which  errors  occur  is  the  high  volume 
expanded  gas  internal  energy  (and  consequently  the  enthalpy),  near  and  above  T  .  For  a 
given  temperature,  the  internal  energy  of  a  single  phase  .as  tends  towarls  that  of  a 
perfect  gas  ns  the  volume  increases.  The  internal  energy  of  a  perfect  gas  i-  only  a 
function  of  temperature  and  consequently  the  high  volume  internal  energy  given  in 
Alderson's  isothermal  tabulation  should  tend  to  the  perfect  gas  limit  nml  not  diver-re  ns 
i-  evident  from  the  tabulation. 

The  effort  of  these  errors  upon  anv  conclusions  of  hypothetical  arc i dent  analysis  is 
riot  expected  to  he  significant. 


APPENDIX  4 


Table 


Table 


Table 


Table 


Table 


Table 


Table 


Tabular  Format  of  the  Equation  of  State 


Saturation  temperature 


listing: 


listing: 


specific  volume 
internal  energies 
enthalpies 
entropies 
pressure 

J 

Saturation  pressures 

' 

specific  volumes 
internal  energies 
enthalpies  ’ 

entropies 
temperatures 

* 

Constant  temperatures 


for  liquid  and  vapour 


for  liquid  and  vapour 


listing:  pressure 

condition 
specific  volume 
internal  energy 
enthalpy 
entropy 


Constant  entropy 


listing:  pressure 

specific  volume 
internal  energy 
enthalpy 
temperature 


Constant  volumes 


listing:  temperature 

internal  energy 
enthalpy 
entropy 
pressure 

Constant  internal  energy 

listing:  temperature 

specific  volume 
enthalpy 
entropy 
pressure 


Saturated  properties 


listing: 


temperature 

thermal  expansion  coefficient,  a,p 
isothermal  compressibility, 
thermal  pressure  coefficient,  yv 
constant  pressure  specific  heat,  C 
constant  volume  specific  heat,  Cv 
Gruneisen  gamma,  I"v 


|  for  liquid  and  vapour 


J 
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APPENDIX  5 


The  Comparison  of  a  Gruneisen  Extrapolation  with  One 
Based  on  the  Thermal  Pressure  Coefficient,  and  the 
Implications  for  the  Behaviour  of  C 

In  contrast  to  the  thermal  pressure  coefficient  teclmique,  a  Gruneisen  extrapolation 
asserts  that  the  Gruneisen  parameter  is  invariant  with  volume.  It  is  defined  as 


Comparing  (5.1)  with  (l.l),  we  see  that  in  using  a  description  we  are  formulating 
an  equation  of  state  in  terms  of  (P,V,E)  rather  than  (P,V,T).  Different  techniques  have 
to  be  employed  to  develop  such  an  equation  of  state;  then,  the  two  extrapolative  tech¬ 
niques  can  be  compared. 

Calculations  for  both  systems  have  been  performed  by  extrapolating  along  isotherms 
away  from  the  saturated  liquid  line  into  the  compressed  liquid.  Fig.  5  shows  isotherms 
for  T  $  2000K.  Below  this  temperature,  in  fact  up  to  ca.  2300K,  the  difference  in  the 
two  techniques  is  very  small  over  the  range  of  pressures  shown  in  Fig.  5.  Beyond  this 
pressure  range  the  isotherms  begin  to  diverge.  Above  2300K  the  isotherms  for  the  two 
techniques  were  clearly  different,  as  can  be  seen. 

The  agreement  between  the  y  and  techniques  is  illustrated  by  the  following:- 
from  (5.1)  we  can  write 


r 

V 


(5.; 


If  yv  and  I"v  are  only  functions  of  V,  it  follows  from  (5.2)  that  Cy  must  also  only  be  a 
function  of  V.  Employing  the  following  relationship 


we  see  that  the  last  term  is  zero  as  y^  is  only  a  function  of  volume;  consequently 


Thus  Cv  must  be  a  constant  if  yy  and  are  functions  of  volume  only.  It  has  been  pre¬ 
viously  demonstrated  (Fig.  2)  that  C  ,  for  both  liquid  and  vapour,  is  reasonably  constant 
up  to  ca.  2-300K.  They  can  now  be  seen  to  be  sufficiently  constant  as  to  introduce  no 
significant  difference  in  the  yy  and  r  techniques  up  to  ca.  2300K.  Above  this  tempera¬ 
ture  C  for  both  phases  rapidly  increases,  and  a  commensurate  divergence  in  y  ,  F 

v  V  V 

extrapolation?  is  observed. 
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Volume  — ► 


A. E  R  E.  R 9847.  FIG.1.  THE  GENERAL  P-V  SURFACE. 
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A.E  R.E.R9847.  FIG.  3.  EVALUATION  OF  THE  INTERNAL  ENERGY  IN 

THE  SUBCRITICAL  COMPRESSED  LIQUID 
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A.E.R.E. R9847.  FIG. 4.  EVALUATION  OF  THE  INTERNAL  ENERGY  IN  THE 

SUPERCRITICAL  COMPRESSED  LIQUID 
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